Astronomy <fc Astrophysics manuscript no. 6845 


© ESO 2008 


February 4, 2008 





A new Neptune-mass planet orbiting HD 219828' 



C. Melo 1 , N. C. Santos 2 ' 3 ' 4 , W. Gieren 5 , G. Pietrzynski 5 , M. T. Ruiz 6 , S. G. Sousa 2 ' 7 ' 8 , 
Lovis 3 , M. Mayor 3 , F. Pepe 3 , D. Queloz 3 , R. da Silva 3 , and S. Udry 3 



F. Bouchy 9 , C. 



o 
o 

(N 

X> 
<D 
PL, 

O 

(N 

(N 
> 

in 

(N 
O 

o 

Oh 
6 



European Southern Observatory, Casilla 19001, Santiago 19, Chile 

Centro de Astronomia e Astroffsica da Universidade de Lisboa, Observatorio Astronomico de Lisboa, Tapada da 
Ajuda, 1349-018 Lisboa, Portugal 

Observatoire de Geneve, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 
Centro de Geofisica de Evora, Rua Romao Ramalho 59, 7002-554 Evora, Portugal 
Universidad de Concepcion, Departamento de Fisica, Casilla 160-C, Concepcion, Chile 
Departamento de Astronomia, Universidad de Chile, Casilla Postal 36D, Santiago, Chile 
Centro de Astroffsica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 
Departamento de Matematica Aplicada, Faculdade de Ciencias da Universidade do Porto, Portugal 
Institut d'Astrophysique de Paris, 98bis Bd. Arago, 75014 Paris, France 



Received XXX; accepted XXX 



ABSTRACT 



Two years ago a new benchmark for the planetary survey was set with the discoveries of three extrasolar planets with 
masses below 2OM0. In particular, the serendipitous discovery of the 14M® planet around fx Ara found with HARPS 
with a semi-amplitude of only 4 ms _1 put in evidence the tremendous potential of HARPS for the search of this class 
of very low-mass planets. Aiming to discovering new worlds similar to /i Ara b, we carried out an intensive campaign 
with HARPS to observe a selected sample of northern stars covering a range of metallicity from about solar to twice 
solar. Two stars in our program were found to present radial velocity variations compatible with the presence of a 
planet-mass companion. The first of these, HD 219828, was found to be orbited by a planet with a minimum mass 
of 19.8 and an orbital period of 3.83 days. It is the 11th Neptune-mass planet found so far orbiting a solar-type 
star. The radial velocity data clearly show the presence of an additional body to the system, likely of planetary mass. 
The second planet orbits HD 102195, has a mass of 0.45Af./ up and an orbital period of 4.11 days. This planet has been 
already announced by Ge et al. (2006). Our data confirm and improve the orbital solution found by these authors. 
We also show that the high residuals of the orbital solution are caused by stellar activity, and use the bisectors of the 
HARPS cross-correlation function to correct the noise introduced by stellar activity. An improved orbital solution is 
obtained after this correction. This kind of analysis may be used in the future to correct the radial- velocities for stellar 
activity induced noise. 



Key words, stars: individual: HD 219828 
formation - techniques: radial-velocity 



1. Introduction 



The serendipitous discovery by ISantos et alj ((2004a) of a 
14M0 planet orbiting /i Ara, simultaneously followed by the 
anno uncement of another pl anet of similar mass around 55 
Cnc (jMcArthur et al.ll2004D and of a 21M e mass planet 
around G J 436 b ( Butler et al.l l2004h set a new bench- 
mark for planet su rveys. Since then, a few more case s have 



been announced dRivera et al 



IVoet et al.ll2005t lUdrv et al 



three Neptune-mass bodies (jLovis et al.l feOOB). For the first 
time in the literature, the minimum masses of a planet 
found by Doppler shift techniques were indicated in Earth- 
masses and not in Jupiter- or Saturn-masses. Even more 
exciting was the fact that based on their orbital character- 
istics and masses, it is strongly s uggested that these planets 
may be rocky or icy in nature (jAlibert et al. I [20061 ). These 



* Based on observations collected at the La Silla Parana 
Observatory, ESO (Chile) with the HARPS spectrograph at 
the 3.6m telescope (Observing runs 075.C-0332, 076.C-0155, and 
77.C-0101). 



stars: individual: HD 102195 - stars: planetary systems - planetary systems: 



discoveries may have thus unveiled the first super-earths, 
and shown that telluric planets may be common in the so- 
lar neighborhood. 

The detection of these kinds of planets is not easy to be 
achieved. For example, the discovery of the 14M ffi compan- 
ion of n Ara was only made possible thanks to the two main 
factors. First, the intrinsic stability of an instrument like 
HARPS which is able to keep a long-term (several years) 
radial velocity accuracy better than 1 ms . Second, the 
observational strategy, /i Ara was observed in a context of a 
asteroseismology program, where a unique target was mon- 
itored during the whole night and a large number of spectra 
was collected with exposure times of about a 2 minutes, in 
order to measure the stellar acoustic oscillation modes. Seen 
from the perspective of a radial velocity planetary search, 
these oscillations are unwanted and are considered as stellar 
noise. This intrinsic noise is well illustrated by the case of 
H Ara: during the nights when the star was measured with 
HARPS over a several hours period, the residuals around 
the orbital fit were merely around 0.4ms -1 , a value that in- 
creased to near 1.5 m s _1 for the remaining nights, when the 



. .20051 : iBonfils et al.l 120051: 
20061). inc luding a system of 
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radial-velocity of the star was an averag e of only about 15 
minute integrations (|Santos et a l. 2004a); this value would 
be around 3 m s _1 if only one single (short exposure) radial- 
velocity measurement was done. In other words, to be able 
to achieve the 1ms -1 precision one needs to average the 
radial velocities for a given star over periods of several min- 
utes (at least 15-20 minutes), making the searches for very 
low mass planets very t ime consuming (see the seismologi- 
cal results on \i Ara by iBouchv et aO2005f) . 

Surprisingly enough, /i Ara and 55 Cnc were back then 
the only two stars that were (incidentally) measured with 
enough precision so that such low mass planets could have 
been detected. This may imply that very low mass short 
period planets are v ery common. Th i s result is support ed by 
recent simulations (jlda k. Linll2004UBenz et alJ l2006) that 
suggest that very low-mass planets may be more frequent 
than the previously found giant worlds. 

These theoretical results and the discovery of several 
Neptune-mass planets suggested that if solar type stars 
are correctly monitored, the number of detection of such 
low-mass rocky planets can be very high. Having this in 
mind, we started a project to use the HARPS spectro- 
graph to monitor a sample of about 40 northern stars with 
metallicity ranging from solar to about twice solar. In the 
present paper, we present the discovery of a new Neptune- 
mass planet orbiting the star HD 219828, as well as the 
confirmation of a Jovian- mass planet orbiting HD 102195. 
We also find evidence for the existence of a longer period 
planet companion to the former star. The presentation of 
the whole sample, along with a detailed study of the be- 
havior of the radial velocities with other properties (e.g. 
spectral type and activity level) is postponed to a future 
paper. 



2. Observations 

The observations were carried out using the HARPS spec- 
trograph (3.6-m ESO telescope, La Silla, Chile), in three 
different observing runs between May 2005 and July 2006 
(ESO observing runs 075.C-0332, 076.C-0155, and 77.C- 
0101). Further measurements of the two stars discussed in 
this paper were also done in different HARPS GTO runs, 
in collaboration with the Geneva team. 

Radial-velocities were o btained directly using the 
HARPS pipeline. We refer to lPepe et ail (|2004h for details 
on the data reduction. The individual spectra were also 
used to derive both the Bisector Inverse Slope (BIS) of the 
HARPS Cross-Corr elation Function (CCF), as defined by 
iQueloz et ail (|200lD . as well as a measurement of the chro- 
mosphe ric activity index log , following the same recipe 
used bv lSantos etal] (|2000f ) for CORALIE spectra. Finally, 
the combined high S/N HARPS spectra were analyzed to 
derive stellar atmospheric para meters and iron abu ndances 
using the method described in lSantos et al.l (|2004bl ). 



Table 1. Stellar parameters for HD 219828. 



3. A new planet around HD 219828 

3.1. Stellar characteristics 



Parameter 


Value 


Reference 


Spectral type 


G0IV 


Hipparcos 


Parallax [mas] 


12.33 ± 1.01 


Hipparcos 


Distance [pc] 


81.1 


Hipparcos 


m v 


8.04 


Hipparcos 


B - V 


0.654 


Hipparcos 


M v 


3.49 


— 


Luminosity [Lq] 


3.34| 




Mass [Mq] 


1.24 


Schaerer et al. (1993) 


log R' HK 


-5.04 


HARPS 


v sini [km s _1 ] 


2.9ft 




T off [K] 


5891 ± 18 




log g 


4.18 ± 0.02 




ft 


1.18 ± 0.02 




[Fe/H] 


+0.19 ± 0.03 





f Using the bolometric correction of IFlowerl l|l996f ) 

|| From HA RPS spectra using a calibration similar to the one 

presented bv lSantos et al.l (|2002T ) 



According to the Hipparcos catalog (|ESAlll997f ). HD 219828 
is a G0IV star with a parallax ir — 12.33 ± 1.01 mas, 
an apparent magnitude m v = 8.04, and a colour index 
B — V = 0.654. From these parameters, and taking the 



bolometric correction from IFlowerl (j 19961 ) . we derive an ab- 
solute magnitude M v =3.49 and a luminosity of 3.34 Lq. 

The analysis of our combined HARPS spectra, with a 
total S/N ratio above 500, provide T eff = 5891 ±18, log g = 
4.19±0.02, and [Fe/H] = 0.19±0.02. Along with the above 
mention ed absolute magnitu de, the Geneva Evolutionary 
models (jSchaerer et al.lll993[ ) indicate that HD 219828 has 
a mass of 1.24M Q and an age of 5.2 Gyr. In Table [T] we 
summarize the stellar parameters derived for this star. 

HD 219828 shows a slightly lower spectroscopic gravity 
as compared to main-sequence dwarf stars of the same tem- 
perature. An even lower gravity (log g = 4.03) is obtained 
if the isochrone mass is used in conju nction with the stel- 
lar parallax fe.g. ISantos et al1l2004bf ). This result is com- 
patible wit h the spectral classification of the star (G0IV - 
lESAlll997h . and indicates that HD 219828 may be slightly 
evolved. The stellar radius estimated using the luminosity- 
temperature-radius relation is 1.76i?0. 

From the HARPS spectra we derived both a chromo- 
spheric activity index (logi?^ K = —5.04) and an estimate 
of the projected rotational velocity of the star (v sini = 2.9 
km s _1 ). From the activity level and the B — V colour 
we derive an age of 6.5 Gyr ([Henry et al.l [l996) (at least 
above 2 Gyr - iPace fc Pasquinill2004D . and a rotational pe- 
riod of 26 days (jNoves et al.l Il984f ). All these values sug- 
gest that HD 219828 is an old chromospherically quiet star. 
According to the Hipparcos catalog, this star is consid- 
ered to be constant in photometry, with a scatter of only 
0.013 mag (typical for a constant star of its magnitude). 

3.2. HARPS orbital solution 

HD 219828 was observed 22 times with the HARPS spectro- 
graph, between May 2005 to August 2006. Each measure- 
ment was done using an exposure time of 90 s, in order to 
avera ge out the stellar oscillation noise (e.g IBouchv et al.l 
I2005D . The complete radial velocity measurements obtained 
and the errors are presented in Tabled It is worth noticing 
that the errors quoted in Tabled which are used to plot 
the error bars, refer solely to the instrumental (calibration) 
and photon-noise error share of the total error budget (e.g. 
activity and/or stellar oscillations are not considered, given 
the difficulty in having a clear estimate of their influence). 
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Table 2. Radial-velocity measurements of HD 219828 and 
HD 102195. 

JD V r [kms' 1 ] a(V r ) [kins' 1 ] = 



HD 219828 








2453509 


928056 


-24 


03256 


0.00066 


2453510 


928367 


-24 


01790 


0.00091 


2453701 


542843 


-24 


02993 


0.00041 


2453702 


532514 


-24 


01817 


0.00046 


2453703 


531084 


-24 


01824 


0.00036 


2453704 


530683 


-24 


02639 


0.00040 


2453705 


544423 


-24 


02791 


0.00035 


2453706 


530181 


-24 


01890 


0.00045 


2453707 


529860 


-24 


02140 


0.00059 


2453708 


551863 


-24 


03220 


0.00046 


2453709 


539910 


-24 


02877 


0.00057 


2453710 


539622 


-24 


01823 


0.00049 


2453930 


824911 


-24 


05650 


0.00066 


2453931 


779782 


-24 


05201 


0.00104 


2453932 


776323 


-24 


04571 


0.00065 


2453933 


735140 


-24 


04473 


0.00080 


2453934 


739643 


-24 


05830 


0.00070 


2453935 


739712 


-24 


05720 


0.00069 


2453936 


730288 


-24 


04521 


0.00076 


2453946 


762931 


-24 


06075 


0.00071 


2453951 


833961 


-24 


05168 


0.00074 


2453975 


734459 


-24 


05942 


0.00069 


HD 102195 








2453501 


574413 


2.15457 


0.00037 


2453503 


580640 


2.06803 


0.00133 


2453504 


587954 


2.15706 


0.00082 


2453506 


621511 


2.08172 


0.00146 


2453550 


535485 


2.17258 


0.00075 


2453551 


523110 


2.11355 


0.00117 


2453757 


789800 


2.06038 


0.00038 


2453761 


827482 


2.05137 


0.00036 


2453765 


771759 


2.07657 


0.00037 


2453785 


772065 


2.13058 


0.00037 


2453788 


796126 


2.16390 


0.00033 


2453927 


478651 


2.08383 


0.00062 


2453930 


461429 


2.05648 


0.00062 


2453931 


470554 


2.06908 


0.00067 


2453932 


485026 


2.17049 


0.00065 


2453933 


465358 


2.18322 


0.00069 


2453934 


470134 


2.08446 


0.00062 


2453935 


471600 


2.06340 


0.00065 


2453936 


476701 


2.13495 


0.00061 



A quick analysis of the data revealed the presence of 
a stable 3.83-day period radial-velocity signal. This signal 
is best fitted using a Keplerian fit with an amplitude K 

Table 3. Elements of the fitted orbit for HD 219828b. 



p 


3.8335±0.0013 


fl 


T 


2453898.6289±0.072 


el 


e 


0.0 f 




V r 


-24.032±0.001 


kms -1 


U) 


0.0 f 


degr 


Ki 


7.0±0.5 


ms -1 


Mm) 


0.1364 ■ 10~ 12 ±0.3572 ■ 10~ 12 


M 


o(0 - C) 


1.7 


ms 1 


N 


22 




J7l2 sini 


19.8 


M ffi 



' Fixed to zero. The obtained eccentricity was consiste nt with a 
circular orbit according to the lLucv fc Sweeney! f| 1 9 7 lh test. 



3 




' i i i i i i i i i i i i i i i i i i i i i i_ 

100.0 102.0 104.0 106.0 108.0 110.0 
JD (-2453600.0) 



Fig. 1. Top: Phase-folded radial- velocity measurements of 
HD219828, and the best Keplerian fit to the data with a 
period of 3.8-days. In this fit the long term radial-velocity 
trend was subtracted. See text for more details. Bottom: 
Radial-velocity measurements of HD219828 as a function 
of time for one epoch of measurements. 



of 7.0ms _1 , and a non-significant eccentricity (see Fig. [2 
leaving the eccentricity as a free parameter in the orbital 
fit provides a value of 0.11±0.09). This corresponds to the 
expected signal induced by the presence of a 19.8 Earth- 
mass (minimum-mass) companion to HD 219828 (Tabled]). 

Our radial velocity measurements as a function of Julian 
Date are shown in FigO The plot shows that, in addition 
to the short period signal, there is a long-period trend su- 
perimposed. We have tried to fit a two Keplerian model 
to the whole radial- velocity set using the stakanof genetic 
algorithm, recen tly developed by Tamuz et al. (in prepa- 
ration; see also iPepe et al]|2006[ ). Although the 3. 83-day 
period signal was always present in the solutions, several 
similar quality solutions were found for the long period 
signal. These always ranged from ^180 to ~800 days in 
period, corresponding to the presence of a planet in the 
Jupiter-mass domain (^0.7Mj up ). The best fit (with an 
rms of only 1.2 ms -1 ) gave a period of 181-days, eccentric- 
ity of 0.3, and amplitude K=21.6ms~ 1 for the long period 
orbit. 

Given the ambiguity in the result, we prefer at this point 
to fit a simple quadratic drift to the residuals of the short 
period signal. The global rms of the Keplerian + drift so- 
lution is 1.7ms _1 . We note, however, that this long term 
trend does not perfectly fit the observed variation; some of 
the two Keplerian models had better x 2 than the adopted 
preliminary solution. This is also illustrated by the increas- 
ing residuals of the fit as a function of time (Fig. [2]). These 
facts explain the relatively high rms of the orbital solution 
found, and may imply that the orbital parameters (e.g. e 
and T) of the short period planet may change slightly when 
the long period solution is better constrained. 
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-100 100 200 300 -10 -5 5 10 

JD (-2453600.0) RV [m/s] 



Fig. 2. Top: Radial-velocity measurements of HD219828 as 
a function of time, and a fit to the data including a 3.8- 
day period Keplerian and a long period quadratic trend. 
Bottom: Residuals of the fit. 



To understand if the short period radial-velocity signal 
observed could be due to the presence of stellar spots (e.g. 
ISaar fc Donahudll997t iQueloz et all200lD or stellar blends 
(|Santos et alj l2002t ). we computed the Bisector Inverse 
Slope (BIS) of the HARPS Cross-Correlation Functions 
(CCF) of HD 219828. In Fig.[3] we plot the resulting BIS 
as a function of the radial- velocity, after having subtracted 
the long period quadratic trend to the data. The result 
shows that no correlation exists between the two quanti- 
ties, suggesting that stellar activity of stellar blends cannot 
explain the short period and low amplitude radial-velocity 
variation observed. 

Together with the very low activity level of the star, 
we thus conclude that the 3.83-day orbital period observed 
can be better explained by the presence of a Neptune-mass 
planet orbiting HD 219828. 



4. The planet around HD 102195 

4.1. Stellar characteristics of HD 102195 



According to the Hipparcos catalog (|ESAlll997f ). HD 102195 
is a KOV star with visual magnitude m v = 8.07 and colour 
index B — V — 0.835. From the visual magnitude and as- 
trometric distance of 28.98 pc (it = 34.51 ± 1.16), we derive 
an absolute magnitude M v = 5.76. Consid ering the bolo- 
metric correction taken from iFlowerl (| 19961 ) . the obtained 
luminosity is then 0.47 L©. 

Accurate stellar atmospheric parameters for HD 102195 
were derived from the analysis of our combined HARPS 
spectra, with a S/N ratio above 300. Our analysis provides 
T eff = 5291 ± 34K, log g = 4.45 ± 0.04, and [Fe/H] = 
0.05 ± 0.05, in e xcellent agreement with previous estimates 
(|Ge et al.ll2006f) . This value for the gravity is compatible 
with the one derived using the temperature, astrometric 



Fig. 3. BIS vs. radial-velocity for HD 219828. The best lin- 
ear fit to the data is shown, and has a slope compatible 
with zero. 



distance, and stellar luminosity (logg = 4.56). A stellar 
mass of 0.87M© is derived from the absolute magnitude, ef- 
fective temperature, and metallicity, after com parison with 
stellar evolution models (|Schaller et al.l 1992V The stellar 
radius estimated using the luminosity-temperature-radius 
relation is 0.82i?©. According to the properties summarized 
in Table |U HD 102195 is a typical solar metallicity main- 
sequence dwarf. 

The HARPS cross-correlation function gives an esti- 
mate of the projected stellar rotational velocity v smi = 
2.6 km s -1 , a va lue that is reas onably compatible with the 
one derived by (|Ge et al .|[2006) . Using HARPS spectra, we 
further derive a rather high value for the chromospheric 
activity index for this star (logi?^ K = —4.56). A sim- 
ilarly high value of lo gi?^ = —4.30 was obtained by 
IStrassmeier et al.1 (|200Clh . 

Using the derived chromosp heric activity index , and the 
colour of the star, we used the lHenrv et al.l (|l996f ) calibra- 
tion to derive an estimate for the stellar age of 1.17 Gyr, 
and a rotational period of about 20-days (|Noves et alj 
Il984t ). slightly above the val ue found using the photomet- 
ric light curve (-12-days - iGe et all 12009) . As an addi- 
tional age indicator, the Li abun dance was computed as de- 
scribed in llsraelian et al.l (|2004|). We found a Li abunda nce 
log e (Li) < 0.2. According to ISestito fc Randichl (|2005l ). a 
T e ff ~ 5300K star with its lithium content should have 
an age above 0. 5 Gyr, also compatible with the estimate of 
IGe et all (120061) . 

4.2. The HARPS orbital solution 

The discovery of a giant planet around HD 102195 has been 
recently announced by the Exoplanet Tracker tean{3 using 
a new concept of instrument combining spectroscopy and 
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Table 4. Stellar parameters for HD 102195. 



Parameter 


Value 


Reference 


Spectral type 


KOV 


Hipparcos 


Parallax [mas] 


34.51 ± 1.16 


Hipparcos 


Distance [pc] 


28.98 


Hipparcos 


m v 


8.07 


Hipparcos 


B-V 


0.835 


Hipparcos 


M v 


5.76 


— 


Luminosity [Lq] 


0.47t 




Mass [M&] 
log R' HK 


0.87 


Schaller et al. C1992) 


-4.56 


HARPS 


v sini [km s" 1 ] 


2.6ft 




T efi [K] 


5291 ± 34 




logs 


4.45 ± 0.04 






0.89 ± 0.05 




[Fe/H] 


0.05 ± 0.05 





f Using the bolometric correction of iFloweil (|1996| ) 

ff From HA RPS spectra using a calibration similar to the one 

presented bv lSantos et alj (|2002T ) 




-100.0 0.0 100.0 200.0 

JD (-2453600.0) 



300.0 



Fig. 4. Top: Phase-folded radial-velocity measurements of 
HD 102 195, and the best Keplerian fit to the data with a 
period of 4.1-days. Bottom: Residuals of the Keplerian fit 
as a function of time. 



interferometry (|Ge et all [2006). Since this star is among 
our surveyed sa mple, it is intere sting to compare the re- 
sults obtained bv lGe et alJ (|2006f ) with those obtained with 
HARPS. 

We observed HD 102195 four times in May 2005, in 
the framework of our program aimed at finding earth-mass 
planets with HARPS around a metallicity biased sample of 
stars. Based on these first four measurements, HD 102195 
was clearly detected to be a radial velocity variable. These 
first four measurements did not show any correlation be- 
tween the bisector and the radial velocities, suggesting that 
this star could be indeed a planet host. 

The HARPS GTO team was contacted and kindly ac- 
cepted to monitor HD 102195. In total, and counting also 



S 




RV [m/s] 



Fig. 5. BIS vs. radial-velocity for HD 102195. The best lin- 
ear fit to the data is shown, and has a slope compatible 
with zero. 




o 10 
Residuals [m/s] 



Fig. 6. BIS as a function of the residuals of the 4.1-day 
Keplerian radial-velocity fit for HD 102195. The data is 
clearly correlated. The best linear fit has a slope of —0.73. 



the data points obtained in our subsequent observing runs, 
19 radial- velocity points were added between June 2005 and 
February 2006 in order to derive a good orbital solution. In 
th e meanwh i le, th e star was announced as a planet host 
by iGe et all (|2006h . In Table [2] we show the radial velocity 
measurements collected for HD 102195 between May 2005 
and February 2006 covering a time span of about 300 days. 
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Table 5. Elements of the fitted orbit for HD 102195b, using the original and the corrected radial-velocity measurements. 
See text for more details. 







V_y(JI ILL LULL VLlULlLlLh 




D 
J 


4.114zzo±U.UUUoOO 




a 


T 


2453896.00±0.04 


2453895.96±0.03 


d 


e 


o.ot 


O.Of 




V r 


2.113±0.002 


2.132±0.001 


kms -1 


U! 


O.Of 


O.Of 


degr 


Kj 


64±3 


63±2 


ms -1 


Mm) 


1.135 ■ 10~ 10 ±0.169 ■ 10~ 10 


1.070 ■ 10- 10 ±0.103 • 10 _1 ° 


M 


a(0 - C) 


9.4 


6.1 


ms 1 


N 


19 


19 




mi sini 


0.46 


0.45 


Mj„ p 



Fixed to zero. The obtained eccentricity was consistent with a circular orbit according to the lLucv fc S weeney (1971]) test. 



Using the radial velocity data shown in Table [2] an or- 
bital solution was found (Table [5]). The best fit Keplerian 
orbit gives a period of 4.11-days, an amplitude K of 
64 ms -1 , and a non-significant ecc entricity. These values 
are similar to the ones derived by iGe et alj {2006). The 
phase-folded radial velocity plot of the measurements is 
shown in Fig. [¥] 

Despite the high activity level of the star, an analysis 
of the bisector of the HARPS cross-correlation functions 
shows no clear correlation with the radial- velocities (Fig.0 . 
This strongly suggests that the radial- velocity signal is not 
being produced by stellar spots or stellar blends, and sup- 
ports the planetary explanation. 

Using the best-fit orbital solution and the mass of 
HD 102195, we derive for the companion a minimum mass 
0.46M,7 MD , in cl ose agreement with the 0.49Mj up found by 
IGe et all (120061, . 

The stellar radius together with the rotational period 
mentioned above provide a value for the stellar rotational 
velocity of ~3.5kms . Comparing this with the esti- 
mated value for the projected rotational velocity (v sini = 
2.6kms~ 1 ) we get s'mi = 0.74, corresponding to an inclina- 
tion of 47 degrees. If the orbital axis of the pl anetary orbit 
is ali gned with the stellar rotation axis (e.g. iQueloz et al.l 
l2000f ). such an inclination coul d explain why no transit 
of HD 102195 b was detected by IGe et al.l (|2006l ). The real 
mass of planet would then be 0.62 Mj up . 

4.3. Correcting the radial velocities for the stellar noise 

In the case of HD 102195, the residuals (O-C) given in 
Table[5]of 9.4 m s -1 are much higher than the errors on ra- 
dial velocity quoted in Tabled This extra- noise is lik e ly due 
to stellar activity (see previous section). IGe et al.l ((2006) 
found that this star has a photometric variation with a pe- 
riod of ~12 days, and an amplitude of ~0.015 magnitudes. 

The signature of its high activity level is also seen on the 
bisector analysis done in this paper. Figure [Fj] shows that 
the residuals of the best Keplerian orbit (previous section) 
do correlate with the bisector measurements in a clear indi- 
cation that the high-activity level o f this star is responsible 
for the observed h igh residuals (e.g. lSaar fe Donahud[T997l 
ISantos et ai1l2000f) . 

According to Fig. the relation between the residual 
radial-velocities and the BIS measurements can be approx- 
imated by a linear function a x BIS + b whose coefficients 
were determined by a least-squares fit. This function can 



be subtracted from the radial velocity data points in order 
to correct for the stellar activity. 

Using those corrected velocities, a new orbital solution 
was found. The results, also presented in Tabled show a 
clear improvement. Not only the rms around the orbital fit 
was reduced from 9.4 to 6.1ms _1 (equivalent to subtract- 
ing quadratically a noise of 7ms _1 ), but also the estimated 
errors in all orbital parameters decreased by almost a fac- 
tor of two. The resulting mass for HD 102195 b also slightly 
decreased (0.45 Mj up ) due to the small decrease in the am- 
plitude of the orbital solution. 

This result shows that the study of the BIS can be used 
to correct, at least to some extent, the noise on the radial- 
velocities introduced by intrinsic stellar activity. If the BIS 
vs. RV correlation can be calibrated (observationally or the- 
oretically) as a function of the stellar properties, we could 
think that in the future an externally determined BIS vs. 
RV function can be used to correct the radial velocity data. 
This will be extremely important in reaching accuracy of 
the other of t he cms -1 as expecte d in future instruments 
like CODEX (|Pascmini et al.ll2005fh 

5. Concluding remarks 

In this paper, we present the discovery of a new 20 M ffi 
(minimum- mass) short period planet, orbiting the sun- like 
star HD 219828. The presence of a second longer period 
companion to the system, likely a Jovian planet, is also 
clear from our data. Unfortunately, the bad phase cover- 
age of the radial-velocity points does not permit to settle 
the orbital parameters of this latter companion. Finally, we 
confirm the existence of a m-i sinj=0.45Mj up planet orbit- 
ing HD 102195, previously announced by another team. 

The planet around HD 219828 is the 11th found with 
masses similar to the mass of Neptune. As the number of 
these systems increa ses, new statistical analysis will be pos- 
sible. In this sense, lUdrv et ail (|2006f ) have discussed the 
possibility that the well known strong correlation between 
the presence of planets and the stellar me tallicity that ex- 
ists for stars host i ng gian t planets (e.g. iGonzalezj 119971 : 
ISantos et all l2001l l2004bt iFischer fc Valentil I2005D does 
not seem to be present for their lower mass counterparts. 
For very-low mass companions, the metallicity distribu- 
tion could rather be flat, something that may be explained 
by recent mode ls of planet formation (|Ida fc Linl [2004; 
iBenz et al.ll2006lh Although the metallicity of HD 219828 
is rather high, we must be careful when analyzing the num.- 
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Fig. 7. Transit depth expected for three different planet 
models as a functi on of the mass . The lower curve repre- 
sents the model of I Valencia et al.l ([2006), and the two up- 
per curves the scaled radii considering the mean density 
of Neptune and the icy Saturn satellite Rhea. The vertical 
dotted line represents the location of HD 219828 b. 



bers, since the planet-search project described in this paper 
concentrated its efforts in studying a sample of metal-rich 
stars only. 

An important question triggered by the previous dis- 
coveries is whether these small planets are constituted by 
rocks, ice or a mixture of both. A more solid assessment of 
the nature of HD 219828 b can only be given if the informa- 
tion about its mass and radius is available. The detection 
of a transiting signature would give us information about 
the radius and mean density of the planet, thus opening 
the possibility to unveil its composition, like is the case 
for the giant planets found to transit their host stars (e.g. 
iKonacki et alll2003t [Pont et al.ll2004 ). 

Although the short period of the orbit implies a rea- 
sonable value for the transit probability (^15%), the mea- 
sure of a transit of HD 2 19828 b is not a n easy task. Using 
the models presented by I Valencia et al.l (|2006l ). considering 
that this 77T-2 sin i = 19.8 M ffi planet could have an Earth- 
like composition, we derive its radius to be ~2.2 Earth radii. 
Although this value must be taken as a lower limit (both 
due to the lower limit of the mass, an d to the fact tha t 
it may have an extended atmosphere - lLovis et al.1 [20061. 
observing the transit of such a low-mass planet is likely 
impossible with current ground-based instrumentation (see 
Fig.UJ. 

For the case of HD 102195, we have also tested, for the 
first time, the use of the bisector for the cross-correlation 
function to correct the radial- velocity measurements for the 
noise introduced by stellar activity. The results show that 
a considerable improvement can be found. Although the 
technique may probably be improved, this conclusion con- 



stitutes a first step towards the subtraction of the radial- 
velocity noise produced by stellar activity phenomena. 
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